The structure of matrices in rational Gauss quadrature
نویسندگان
چکیده
This paper is concerned with the approximation of matrix functionals defined by a large, sparse or structured, symmetric definite matrix. These functionals are Stieltjes integrals with a measure supported on a compact real interval. Rational Gauss quadrature rules that are designed to exactly integrate Laurent polynomials with a fixed pole in the vicinity of the support of the measure may yield better approximations of these functionals than standard Gauss quadrature rules with the same number of nodes. It therefore can be attractive to approximate matrix functionals by these rational Gauss rules. We describe the structure of the matrices associated with these quadrature rules, derive remainder terms, and discuss computational aspects. Also, rational Gauss-Radau rules and the applicability of pairs of rational Gauss and Gauss-Radau rules to computing lower and upper bounds for certain matrix functionals are discussed.
منابع مشابه
Rational Gauss Quadrature
The existence of (standard) Gauss quadrature rules with respect to a nonnegative measure dμ with support on the real axis easily can be shown with the aid of orthogonal polynomials with respect to this measure. Efficient algorithms for computing the nodes and weights of an n-point Gauss rule use the n × n symmetric tridiagonal matrix determined by the recursion coefficients for the first n orth...
متن کاملMatrices, Moments, and Rational Quadrature
Abstract. Many problems in science and engineering require the evaluation of functionals of the form Fu(A) = uT f(A)u, where A is a large symmetric matrix, u a vector, and f a nonlinear function. A popular and fairly inexpensive approach to determining upper and lower bounds for such functionals is based on first carrying out a few steps of the Lanczos procedure applied to A with initial vector...
متن کاملComputing rational Gauss-Chebyshev quadrature formulas with complex poles
We provide a fast algorithm to compute arbitrarily many nodes and weights for rational Gauss-Chebyshev quadrature formulas integrating exactly in spaces of rational functions with arbitrary complex poles outside [−1, 1]. This algorithm is based on the derivation of explicit expressions for the Chebyshev (para-)orthogonal rational functions.
متن کاملQuadrature rules for rational functions
It is shown how recent ideas on rational Gauss-type quadrature rules can be extended to Gauss-Kronrod, Gauss-Turr an, and Cauchy principal value quadrature rules. Numerical examples illustrate the advantages in accuracy thus achievable. 0. Introduction The idea of constructing quadrature rules that are exact for rational functions with prescribed poles, rather than for polynomials, has received...
متن کاملComputing rational Gauss-Chebyshev quadrature formulas with complex poles: The algorithm
We provide an algorithm to compute arbitrarily many nodes and weights for rational Gauss-Chebyshev quadrature formulas integrating exactly in spaces of rational functions with complex poles outside [−1, 1]. Contrary to existing rational quadrature formulas, the computational effort is very low, even for extremely high degrees, and under certain conditions on the poles it can be shown that the c...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید
ثبت ناماگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید
ورودعنوان ژورنال:
- Math. Comput.
دوره 82 شماره
صفحات -
تاریخ انتشار 2013